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We study the thermalization process of charm quarks in hot and dense matter. The diffusion of 
heavy quarks is calculated via a Langevin equation, both for a static medium as well as for a QGP 
medium generated by a (3+l)D hydrodynamic model. We define two criteria for the thermalization 
of the heavy quarks, and observe thermalization times that are longer than the lifetime of the QGP 
phase for reasonable values of the diffusion constant. 
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I. INTRODUCTION 

The first decade of RHIC operations have yielded 
a vast amount of interesting and sometimes surprising 
data. It is now generally accepted that RHIC has cre- 
ated a hot and dense state of deconfined quark-gluon 
matter with properties similar to that of an ideal fluid 
- this state of matter has been termed the strongly in- 
teracting Quark- Gluon- Plasma (sQGP) Among the 
unexpected observations to contribute to the moniker 
"strongly interacting" are the surprisingly large values 
of the elliptic flow coefficient V2 and the surprisingly 
small values of nuclear modification factor Raa @t11 ex- 
hibited by D mesons, which are thought to be due to 
charm quarks interacting strongly with the QGP prior 
to hadronization [l^ . 

In principle, the observed behavior of D mesons could 
be indicative of charm quarks thermalizing in the QGP. 
It is therefore of great interest to study the dynamics of 
heavy quarks in an expanding QGP and to verify whether 
they indeed approach equilibrium on a timescale on the 
order of the lifetime of the QGP. Due to the "dead-cone 
effect" [nl, [lH, heavy quark interactions with the QGP 
are dominated by quasi-elastic scattering off light quarks 
and gluons [H, H^- Radiative contributions to e.g. the 
energy-loss exist, but have been shown to be small (10- 
20%) for heavy quarks at small and intermediate trans- 
verse momenta |15| . In the limit of multiple interactions, 
the heavy quark propagation through a QGP can there- 
fore be treated in the framework of a Langevin equatiorQ 
d, IT8l - [2H . In [§] it was estimated that the charm quark 
thermalization time is of the order of 7 fm/c. It was 
also suggested that the presence of resonant heavy-light 
quark interactions at moderate QGP temperatures accel- 
erates the kinetic equilibration of charm quarks [l^ [22| 
and reduces the thermalization time from over 30 fm/c 
(as estimated by previous pQCD calculations [H,[2^) to 
only a few fm/c. 

In this paper we shall use a Langevin equation to 
conduct an investigation into heavy quark thermaliza- 
tion for infinite QGP matter at fixed temperature and 



^ Note that the Langevin equation may not be a good approxima- 
tion if the momentum transfer during each interaction is large. 
See [9lll6l. [T7ll for a discussion on this topic and other constraints 
on the Langevin approach. 



then compare our findings to a dynamical scenario in 
which the heavy quarks propagate through an expand- 
ing and cooling QGP medium, modeled with a (3-|-l)D 
relativistic hydrodynamic approach [23 |. Our goal is 
not to address current experimental data, but to answer 
the question whether general features seen in the data, 
such as the presence of elliptic flow and a small value 
of the nuclear modification factor, can be used to con- 
clude that heavy quarks actually have thermalized in the 
QGP medium created in ultra-relativistic heavy-ion colli- 
sions or whether the heavy quarks remain off-equilibrium 
during their entire evolution, despite exhibiting a strong 
response to the surrounding QGP. 

This paper will be organized as follows. In the next 
section, we will discuss the methodology used in our in- 
vestigation by briefly summarizing the Langevin algo- 
rithm and introducing the criteria we use to test for heavy 
quark thermalization. In Sec. Ill, we will study the ther- 
malization of charm quarks in a static QGP medium at 
flxed temperature and examine how the medium temper- 
ature, initial momentum of the charm quarks and the dif- 
fusion coefficient influence the thermalization process. In 
Sec. IV, the thermalization process in an expanding and 
cooling QGP medium will be simulated. A summary and 
outlook will be provided at the end. 



II. METHODOLOGY 

A. Heavy-quark dynamics 

Due to the "dead-cone effect" , heavy quark interac- 
tions with the QGP are dominated by quasi-elastic scat- 
tering off light quarks and gluons. In the limit of mul- 
tiple interactions, the heavy quark propagation through 
a QGP can therefore be treated in the framework of a 
Langevin equation : 
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In principal, the noise term ^ may depend on the momen- 
tum of the particle, but we assume for simplicity that this 
is not the case here. The random momentum kicks satisfy 
the following correlation relation: 
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In order to simulate the momentum evolution, we 
adopt the Ito discretization: 



p{t + At) = p{t) - dito{m)^t + cAt, 



where 



dito{p) = Vd{p)p- 
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If the transfer of energy is small, it can be shown that the 
fluctuation-dissipation relation applies, which indicates: 
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Furthermore, according to Eq.Q, Gaussian noise with 
width r = k/ At will be used to generate the random 
momentum kicks in our calculation. The diffusion coef- 
ficient is related to the drag term via: 
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B. QGP medium 

For our static infinite matter calculations, the only in- 
formation on the medium that is required is the tempera- 
ture, which remains fixed throughout the time-evolution 
of the heavy quarks. For an expanding QGP, we utilize a 
fully (3-(-l)D rclativistic ideal hydrodynamic model [24 1. 
The initial conditions of the hydrodynamic calculation 
are tuned to describe the hadronic data in the soft sec- 
tor, such as hadron yields, spectra, rapidity-distributions 
as well as radial and elliptic flow. The hydrodynamic 
model provides us with the time-evolution of the spatial 
distribution of temperature and collective flow velocity. 
Using our knowledge of the local flow velocity, for every 
Langevin time step we boost the heavy quark to the local 
rest frame of the fluid cell through which it propagates. 
In the rest-frame of the cell we perform the Langevin 
evolution (at the local temperature of the cell) before 
boosting back to the global CM-frame. 

Among the key assumptions of the hydrodynamic cal- 
culation is that the QCD medium experiences a sudden 
thermalization (to form a QGP) at an initial time tq 
(chosen here to be 0.6 fm/c) at which the hydrodynamic 
evolution commences. The pre-equilibrium phase of the 
medium cannot be modeled via hydrodynamics - here 
we assume that our heavy quarks stream freely prior to 
QGP formation. Any heavy quarks leaving the QGP will 
stream freely as well. In our analysis we only treat de- 
conflncd degrees of freedom and restrict our analysis to 
questions which can be addressed in that context. 



C. Thermalization criterion 

Our thermalization criterion is based on the heavy 
quark energy and its momentum components. For a 
medium at fixed temperature without any inherent col- 
lective flow, an ensemble of thermalized heavy quarks 
has the following energy distribution that allows for the 
straightforward extraction of its "temperature" via an 
exponential fit: 



dN 
pEdE 



(8) 



While this particular form for the energy distribution 
provides a convenient representation for the extraction 
of the "temperature" of the heavy quark ensemble, it is 
insufficient to actually indicate thermalization, since we 
still need to verify isotropy of particle momenta: 

/(k) - C ■ T{^pj + m^ + T)e-^Pi+^-/^. (9) 

Note that if e.g. we initialize our heavy quark ensemble 
with a finite momentum along a given coordinate axis, 
its momentum distribution along that axis will be blue- 
shifted - this can be taken into account by shifting the 
distribution along that axis using a parameter pi] e.g. for 
the z axis this would givc0: 



f{p,) = C ■ T{^[p, - p~,f + to2 + r)e-V(P-P^)=+™VT^ 

(10) 

For an expanding and cooling QGP medium, estab- 
lishing thermalization requires the following procedure: 
for a given time step we select all cells in our hydrody- 
namic medium within a temperature band of T ± AT. 
We then boost all charm quarks localized in those cells 
into the respective local rest frames of the cells they are 
residing in and calculate the resulting heavy quark en- 
ergy and momentum distributions. If the energy and 
momentum distributions of the heavy quark ensemble 
yield thermal distributions with a temperature that lies 
within our selected temperature band, we can conclude 
that the selected heavy quark ensemble has thermalized 
in the medium at the given temperature and selected time 
step. This procedure can then be repeated for other tem- 
perature bands and time steps. 



III. EQUILIBRATION IN A STATIC MEDIUM 

A. Energy vs. momentum distributions as 
thermalization measures 

We start our investigation by simulating the diffusion 
of charm quarks in an infinite medium with a static tem- 



^ Rigorously, should be boosted via 'ypz + 7/3 -E. However, it is 
found that in our study, for /3 not too large (below 0.8), the much 
more convenient Eg. lllOII already fits the spectrum well (with an 
error less than 5% for T). 
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FIG. 1. (Color online) The evolution of the energy spectrum 
with respect to time, (a) shows the results between 2 fm/c 
and 8 fm/c, where no linear relation is observed; and (b) 
shows the results between 10 fm/c and 30 fm/c, where the 
linear relation is apparent. 



pcraturc of 300 MeV. The initial momentum of the charm 
quarks is chosen to be 5 GeV/c in the z direction. The 
diffusion coefficient is set to be Z? = 6/ (27rT), which is the 
value that was found in Q to provide the best agreement 
to data on the elliptic flow and the nuclear modification 
factor of non-photonic electrons [sj . It is worth mention- 
ing that different models for the medium may influence 
the value of the tr ansp ort coefficients required to fit the 
experimental data [2ll |. 

FiglT] shows the energy spectrum dN /pEdE vs. E for 
different diffusion times. We see that a linear relation be- 
tween \n{dN/pEdE) and E docs not occur for diffusion 
times shorter than approx. 10 fm/c. The distribution 
appears thermal for longer diffusion times. However, the 
slope continues to increase as a function of diffusion time 
and does not converge to the temperature of the medium 
until a diffusion time of around 30 fm/c. Therefore, de- 
spite the shape of the energy distribution, our ensemble of 
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FIG. 2. (Color online) The evolution of the momentum spec- 
trum with respect to time. 



charm quarks is not fully equilibrated for diffusion times 
between 10 fm/c and 30 fm/c. 

Figl2] shows the momentum distributions in the three 
directions for the same diffusion times as depicted in 
FiglT] Since the initial momentum of the charm quarks 
is in the z direction, the Px and Py spectra are symmet- 
ric with respect to and are virtually identical to each 
other. The center of the Pz spectrum, which initially is 
a delta function at = 5 GeV/c due to our initial con- 
dition, shifts towards as a function of increasing diffu- 
sion time, signifying the influence of the drag term of the 
Langevin equation on the dominant direction of propa- 
gation. The widths of the momentum distributions along 
all three coordinate axes start to agree with each other 
for diffusion times later than approx. 20 fm/c, hinting at 
a common "temperature" . At that time, full isotropy of 
the momenta is obtained. 

We can obtain the "temperature" evolution of the 
charm quark ensemble by fltting the momentum distri- 
butions with Eq. (|10p . and compare the respective values 
of the temperature parameter with those obtained from 
the slope of the energy distribution. The results can be 
seen in Fig|3]as a function of diffusion time, which helps 
to summarize our main observations from the previous 
figures as follows: 

• for diffusion times less than 10 fm/c, the Pz spec- 
trum is distinctly separated from the Px and Py 
spectra. The temperature parameters extracted 
from the widths of the distributions initially rise 
and are significantly above the actual temperature 
of the medium. They are of different values for the 
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FIG. 3. (Color online) A comparison of "temperatures" ob- 
tained from px, py, and the energy spectra. Since x and 
y directions are symmetric, we take the mean value of the 
"temperatures" fitted from px and py spectra here. 



transverse momentum distributions vs. the longi- 
tudinal momentum distribution (as defined by the 
z axis). No temperature can be obtained from the 
energy distribution since no linear relation is ob- 
served at those short diffusion times. 

• for diffusion times between 10 fm/c and 20 fm/c. all 
momentum distributions as well as the energy dis- 
tribution exhibit a thermal shape, even though the 
extracted temperature parameters strongly differ 
among each other and from the temperature of the 
medium. Interestingly, the temperature parameter 
extracted from the longitudinal momentum distri- 
bution seems to track that from the energy distri- 
bution during this time interval. However, both 
are significantly higher than the temperature pa- 
rameter extracted from the transverse momentum 
distributions. 

• at a diffusion time of roughly 25 fm/c, all temper- 
ature parameters agree well with each other and 
have converged to the temperature of the medium, 
signaling full equilibration of our charm quark en- 
semble. 

Based on our observations, we shall define a "quasi- 
equilibrium" stage to be a near equilibrium state where 
the linear relation of the \n{dN/pEdE) vs. E distribu- 
tion appears and the temperature parameters extracted 
from energy and spectra are of approximately the 
same value. During this stage, our ensemble of charm 
quarks exhibits thermal properties, even though it has 
not yet fully equilibrated with the surrounding medium. 
In contrast, a "full equilibrium" is approached when the 
temperature parameters extracted from different ways 
agree with each other and reach the temperature of the 
medium. 



B. The Blue Shift Effect 

To further understand the nature of the quasi- 
equilibrium state and why the "temperature" of the 
charm quarks is still higher than that of the medium in 
this stage of the evolution, despite the linear relation be- 
tween the \n{dN/pEdE) vs. E distribution, we may have 
to take into account a "blue shift" caused by the center 
of mass motion of the charm quarks: for our analysis, 
all charm quarks were initialized to carry a momentum 
of 5 GeV/c in the z direction. This initial momentum 
can be seen as a center of mass motion - it contributes 
to the total energy of the particle, but is non-thermal in 
origin. As the charm quark diffuses through the medium, 
this initial kinetic energy will dissipate through the drag 
term of the Langevin equation. During this dissipation 
dominated phase of the charm quark evolution, the cen- 
ter of mass motion will contribute to an additional part 
of energy and therefore a higher "temperature" . In this 
sense, the quasi-equilibrium state can be understood as a 
stage when the thermal part of heavy quark motion is al- 
ready close to equilibrium but the center of mass motion 
has not entirely dissipated. 

This blue shift can be verified by fitting Eq. ([TO|) with 
a momentum distribution that includes a momentum 
shift pz. We will see that the blue shift is suppressed 
if the initial momentum of the heavy quarks is small 
enough that it becomes comparable to its thermal mo- 
tion. This in fact helps prove that the blue shift is the 
reason for a higher heavy quark "temperature" in the 
quasi-equilibrium state. 



C. Parameter Dependencies 

Let us now investigate the temperature dependence of 
the thermalization time: as in the previous sections, we 
initialize charm quarks with 5 GeV/c momentum along 
the z direction. The diffusion coefficient is set to be I? = 
6/(27rT). Here, we vary the temperature of the static 
medium and examine its corresponding influence on the 
thermalization process. The results are shown in FigQ] 
and Fig [5] 

In FiglU we compare the thermalization of charm 
quarks in a 300 MeV temperature medium to that in 
a medium at 160 MeV. For charm quarks with an ini- 
tial momentum of 5 GeV/c, their "temperature" ex- 
tracted from the dN/dpz distribution increases first, and 
then decreases until it approaches the medium temper- 
ature at which time the full equilibrium is reached. It 
is observed that the turning point between the rise and 
fall of the "temperature" is close to the onset of the 
quasi-equilibrium stage. A similar trend can be observed 
via the dN/dpx,y distribution (not shown here), but the 
range of variation in the extracted "temperature" is not 
as large as that from the dN/dpz distribution. The above 
observation can also be obtained from media with other 
temperatures indicated in Figl5](not shown here). 
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FIG. 4. (Color online) A comparison of the charm quark piQ q (Color online) A comparison of the thermalization 
thermalization processes in static media with different tem- processes of charm quarks with different magnitude of initial 
peratures. momenta in the z direction. 
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FIG. 5. (Color online) The variation of critical times of ther- 
malization with respect to the medium temperature. 



We can understand the "turning point" behavior as 
follows. Since the charm quarks are all sampled with 
5 GeV initial momentum in the z direction here, they 
are initially at zero temperature with respect to their 
center of mass frame. During the diffusion process, the 
momenta of the charm quarks are randomized and thus 
their "temperature" increases. When the random mo- 
tion approaches equilibrium, i.e., the entrance of the 
quasi-equilibrium state, the charm quarks' "tempera- 
ture" might be higher than the medium temperature if 
the initial center of mass momentum of the ensemble has 
not entirely been dissipated. After that, this tempera- 
ture parameter gradually decreases towards the medium 
temperature with the dissipation of the center of mass 
momentum. Such a rise in the temperature parameter 
before reaching the quasi-equilibrium results from this 
particular initialization of charm quarks, and will not oc- 
cur for more realistic scenarios where charm quarks arc 



initialized in more realistic waysH 

As the temperature of the medium decreases, both the 
time needed to reach quasi-equilibrium and to obtain 
full equilibrium increase. The time to enter the quasi- 
equilibrium stage can vary from 8 fm/c (for 350 MeV 
medium) to 35 fm/c (for 160 MeV medium), and the 
time to approach the full equilibrium can vary from 
20 fm/c (for 350 MeV medium) to 70 fm/c (for 160 MeV 
medium). This variation is clearly shown in Fig IS] Note 
that since the calculation is carried out for static medium, 
critical behavior near Tc ~ 160 MeV is not expected here. 

In order to compare the results for the static medium 
with those for the QGP medium later, we need to sim- 
ulate the diffusion of charm quarks with different initial 
momenta. Due to the expansion of the QGP medium it- 
self, the initial momentum of the charm quarks in the lo- 
cal rest frame might be much smaller than that in the lab 
frame. Therefore, in FiglB]wc show our analysis for charm 
quarks with initial momenta of 1 GeV/c and 3 GeV/c in 
a medium with 160 MeV temperature and compare to 
our previous results. 

We observe that the results of 3 GeV/c initial momen- 
tum are similar to those of the 5 GeV/c case shown in 
FiglD the "temperature" obtained from the dN/dpz dis- 
tribution first increases and then decreases until the full 
equilibrium is approached, and the turning point corre- 
sponds to the entrance of the quasi-equilibrium stage. 
However, the results of the 1 GcV/c initial momentum 



^ The temperature parameter before quasi-equilibrium is ap- 
proached provides some insight into the dynamics of thermaliza- 
tion, namely the interplay of momentum broadening vs. energy 
dissipation through coUisional energy loss. For a fixed initial 
momentum, this interplay is more obvious than compared to an 
initialization with charm quarks of different initial momenta in 
the later sections of our work. However, the values of the "tem- 
perature" obtained in this domain are not indicative of thermal 
behavior but just a measure of the broadening of the charm quark 
momentum distribution, since Eg. lllDI I is only strictly valid near 
equilibrium. 
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FIG. 7. (Color online) The thermalization processes of charm 
quarks with initial momentum in random direction. Since x, 
y and z directions are symmetric in this situation, we take 
the average of the "temperatures" obtained from p^, Py and 
spectra. 
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FIG. 8. (Color online) The variation of critical times of 
thermalization with respect to the diffusion coefflcient of the 
medium. 



case are different. There is no overshoot in the "tempera- 
ture" for this case. Instead, the "temperature" obtained 
from the dN/ dpz distribution keeps increasing gradually 
towards the temperature of the medium until the full 
equilibrium is reached. In other words, unlike in the 
higher initial momentum cases, the "temperature" of the 
charm quarks with 1 GcV/c initial momentum is always 
below that of the medium until reaching full equilibrium. 
We attribute this difference to the suppression of the blue 
shift for the center of mass motion of the charm quarks, 
which is too small in this case to contribute significantly 
to the energy of the charm quark. 

For this low initial momentum situation, the entrance 
to the quasi-equilibrium stage can no longer be deter- 
mined by the turning point, but only by the appearance 
of the linear relation in the \n{dN/pEdE) vs. E distribu- 
tion. In our simulation, for heavy quarks with 1 GcV/c 
initial momentum, this linear relation docs not appear 
until 10 fm/c for the 160 MeV medium. 

Additionally, a comparison between FigU] and FiglH] 
suggests that with a decrease of the magnitude of the 
initial momentum, the thermalization occurs faster. For 
instance, the times needed for 5 GeV/c charm quarks to 
reach quasi-equilibrium and full equilibrium are 35 fm/c 
and 70 fm/c respectively, while these critical times for 
1 GeV/c charm quarks arc 10 fm/c and 20 fm/c. This 
decrease can also be understood with the concept of the 
blue shift: it takes less time for smaller initial momentum 
to dissipate. 

In order to verify the spatial invariancc of our cal- 
culation, we randomize the direction of the initial mo- 
mentum of the charm quarks and examine their ther- 
malization process. The results are shown in Fig|7] for a 
medium temperature of 160 MeV. Comparing with FiglHl 
we observe that even though the trend of the evolution 
of the charm quark "temperature" may appear different, 
due to the coordinate axes now not reflecting anymore 
purely longitudinal or transverse motion respectively, the 



critical times for the charm quarks to enter the quasi- 
equilibrium state and the full equilibrium state are al- 
most identical. For instance, both FiglB] and FiglT] indi- 
cate that for the 160 MeV static medium, charm quarks 
with 3 GeV/c initial longitudinal momentum enter the 
quasi-equilibrium stage at around 25 fm/c, and are en- 
tirely equilibrated at around 50 fm/c. Note that for an 
expanding QGP medium, the initial orientation of the 
motion of the charm quarks with respect the beam axis 
(i.e. the direction of longitudinal expansion) will play a 
significant role. 

To examine the influence of the diffusion coefficient on 
the process of heavy quark thermalization, we initialize 
the medium at a temperature of 300 MeV and the charm 
quarks with an initial momentum of 5 GeV/c in the z 
direction, and vary the diffusion coefficient in accordance 
with the choices in Q. The results arc shown in FiglHl 

The results indicate that as the diffusion coefficient 
decreases, or the drag coefficient increases, the ther- 
malization process can be speeded up obviously. The 
times needed to enter the quasi-equilibrium state and the 
full equilibrium state can be respectively reduced from 
70fm/c and 150fm/c for D = 12/{2ttT) to 10 fm/c and 
25 fm/c for D = 1.5/(27rT). 



IV. ANALYSIS IN A REALISTIC QGP MEDIUM 

In this section, we study the charm quark thermaliza- 
tion process in a realistic expanding QGP medium. The 
QGP medium is generated by a 3-1-1 dimensional ideal hy- 
drodynamic calculation with initial conditions that have 
been adjusted to reproduce bulk properties of the QCD 
medium created in central Au+Au collisions at RHIC 
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FIG. 9. (Color online) The charm quark thermalization pro- 
cess in the QGP medium. The charm quarks have initial 
momentum of 5 GeV/c in the z direction. 



[23|FI If not specified, the diffusion coefficient is chosen to 
he D = 6/(27rr) in this section following the choice most 
compatible with the data in Q. Although after 8 fm/c 
the medium temperature is below Tc (160 MeV), we ex- 
tend our study of charm quarks to the time of freeze-out 
(around 20 fm/c), assuming that D mesons interact with 
the medium in a similar way as the charm quarks do with 
QGP. 

Let us first investigate the evolution of charm quarks 
with 5 GeV/c initial momentum in the z direction. The 
results are shown in FiglHl As can be seen, we are able to 
obtain the linear relation from the \n{dN/pEdE) vs. E 
distribution after 8 fm/c. Meanwhile, the "temperature" 
obtained from the dN/dp^ distribution is always below 
the local temperature of the QGP medium probed by the 
selected sample of charm quarks. 

By examining the value of Pz in Eg. ljlOp . we find that 
the initial momentum of the charm quarks in the local 
rest frame is around 1.8 GeV/c. By comparing with the 
results of charm quarks with 1 GeV/c initial momentum 
in the 160 MeV static medium (shown in FiglHI), we can 
conclude that the charm quarks here enter the quasi- 
equilibrium state after about 8 fm/c, but never reach the 
full equilibrium prior to freeze-out. In other words, the 
charm quarks are approaching the local temperature of 
the QGP medium but never catch up with it. 

To bring the simulation closer to realistic conditions in 
relativistic heavy-ion collisions, where charm quarks are 
distributed mostly throughout the mid-rapidity region, 
we set the initial momentum to be in the z direction, 



Recent study [25j indicates the viscosity of the hydro medium 
may decrease the eUiptic flow of 7r+ and p, but does not have 
a significant influence on their spectra. Likewise, the temporal 
evolution of a viscous hydro (in terms of time-scale, flow veloci- 
ties and temperatures) may differ at most by 20% from the ideal 
hydrodynamic evolution used here - this is well within the sys- 
tematic parametric uncertainties of our study. 
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FIG. 10. (Color online) The charm quark thermalization pro- 
cess in the QGP medium. The charm quarks have initial mo- 
mentum of 5 GeV/c randomly distributed in the transverse 
plane. Since x and y directions are symmetric here, we take 
the mean value of the "temperatures" fitted from px and py 
spectra. 



and randomly distribute its direction in the transverse 
plane with a magnitude of 5 GeV/c. The results are 
shown in Fig llOl Unlike the static medium, the random 
distribution of the initial momentum direction increases 
the complexity of the thermalization process in the QGP 
medium. According to the occurrence of the linear rela- 
tion of the \n{dN/pEdE) vs. E distribution, the critical 
time to enter the quasi-equilibrium state is delayed to 
16 fm/c and the charm quarks remain far off equilibrium 
at the time of freeze-out. 

If we want to see whether thermalization can be 
achieved in principle in an expanding QGP medium, we 
need to increase the interaction strength of the charm 
quarks. This can be realized by decreasing the diffusion 
coefficient to = 1.5/(27rT). For this value, we observe a 
significantly faster thermalization. The results are shown 
in Fig llli which shows that the charm quarks enter the 
quasi-equilibrium stage as early as around 4 fm/c, and 
appear very close to full equilibrium at the time of freeze- 
out. 

Finally, we perform our analysis with charm quarks 
whose initial distribution is sampled from Au-|-Au colli- 
sions at RHIC ener gies generated by the VNI/BMS par- 
ton cascade model |26l . |27[ . The medium evolution via 
a (3-|-l)D hydro and charm quark dynamics remain the 
same as in the previous calculations. The cm. energy 
is ^/s = 200 GeV and the impact parameter is fixed at 
2.4 fm for the generation of both the QGP medium and 
the initial charm quarks. The results are shown in Fig ll2l 



As can be seen in Fig 12(a) with a diffusion coefficient 
of _D — 6/(27rr), the "temperature" of the charm quarks 
never manages to catch up with that of the medium until 
freeze-out. A closer observation indicates that the charm 
quarks remain far off equilibrium during the entire life- 
time of the QGP phase, i.e., when the medium temper- 
ature is above Tc. However, Fig 12(b) suggests that if 
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FIG. 11. (Color online) The charm quark thermalization 
process in the QGP medium with small dilfusion coefficient 
D = 1.5/(27rT). The charm quarks have initial momentum of 
5 GeV/c randomly distributed in the transverse plane. 



the diffusion coefBcicnt is reduced to D = 1.5/{2ttT), 
the thermalization process accelerates significantly and 
the "temperature" of the charm quarks is able to catch 
up with that of the medium during its QGP phase, i.e., 
above Tc- Furthermore, a comparison between Fig |12l 
and Fig llOlfTT] indicates that charm quarks with an ini- 
tial momentum distribution given by the cascade model 
thcrmalize faster than charm quarks with an initial mo- 
mentum of fixed magnitude of 5 GeV/c. This is due 
to the realistic PCM sample containing mostly charm 
quarks in the low momentum region, which are already 
rather close to the thermal momentum scale. 



SUMMARY AND OUTLOOK 
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FIG. 12. (Color online) Thermalization of charm quarks with 
an initial distribution generated by the parton cascade model, 
(a) shows the results in the QGP medium with a diffusion 
coefficient of 73 = 6/{2nT); and (b) shows the results with a 
diffusion coefficient of D = 1.5/(27rr). 



In this paper, we have studied the dynamics of heavy 
quark thermalization in the framework of a Langevin 
equation. Simulations have been carried out in both 
static and dynamic QGP media, and our methodology 
allows for the extraction of the "temperature" of the 
heavy quarks by fitting either the energy or the momen- 
tum spectrum. 

Using an idealized static QGP medium, it is found 
that charm quark thermalization occurs in two distinct 
steps: first a quasi-equilibrium is obtained in which the 
charm quark energy distribution matches that of a ther- 
mal medium, but the momentum distribution remains 
non-isotropic; subsequently the charm quark momenta 
isotropize and the charm quarks are in full equilibrium 
with the surrounding medium. The occurrence of this 
two step process might be explained by the blue shift 
effect due to the center of mass motion of the heavy 
quarks. We define full equilibrium to imply that the 
"temperature" extracted from both methods of fitting 
the energy and the momentum distribution matches that 
of the medium. 



Our simulations in the static medium indicates that as 
the medium temperature decreases, it takes longer time 
for the charm quarks to thermalize. Additionally, the 
thermalization speed is extremely sensitive to the diffu- 
sion coefficient of the medium: as the diffusion coefficient 
decreases, or the drag coefficient increases, the thermal- 
ization process can be significantly accelerated. A de- 
crease in the magnitude of the initial momentum leads 
to a suppression of the blue shift effect and therefore re- 
sults in a faster thermalization process, due to the initial 
charm quark momenta being closer to the thermal mo- 
mentum scale of the system. 

For a realistic expanding QGP medium, we find that 
for choices of the diffusion coefficient that describe the 
elliptic fiow and the nuclear modification factor of the 
charm quarks, thermalization does not occur within the 
lifetime of the QGP phase. Only for a diffusion coefficient 
roughly 1/4 of the desired value do the charm quarks 
actually equilibrate. We find that the manifestation of 
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collective behaviors, such as a significant elliptic flow, or 
the presence of a strongly interacting system (via the nu- 
clear modification factor) is insufficient to conclude that 
charm quarks have actually thermalized in the medium. 
Our investigation indicates that their evolution occurs for 
the most part out of equilibrium, even though their prop- 
erties arc strongly affected by the surrounding medium. 
It is worth emphasizing that this result docs not imply a 
small value as D = 1.5/ (ZttT) for the diffusion coefficient, 
but only indicates that with a currently favored value of 
D = 6/(27rT), which describes experimental results well, 
charm quarks may remain off-equilibrium during their 
entire evolution although they exhibit a strong response 
to the surrounding QGP. 

The various parameter studies undertaken in our anal- 
ysis should provide a measure for the systematic un- 
certainties arising from the initial distribution of charm 
quarks and the variation of the diffusion coefficient. We 
conclude that the charm quarks' deviation from equilib- 
rium during QGP's lifetime is significantly larger than 
one could overcome with different choices for the initial 
state in combination with realistic values for the diffusion 
coefficient around D = 6/{2ttT). Of course, the micro- 
scopic structure of the charm quark scattering cross sec- 
tion in the non-perturbative regime remains largely un- 
known - for example, it has been shown in [Tol . that 
the introduction of resonant heavy-light quark interac- 
tions is able to speed up the heavy quark thermalization 



process. It has also been suggested that the memory ef- 
fect due to strong correlations in the medium [28l - l30j may 
accelerate the thermalization process. In addition, most 
current calculations do not account for any interaction in 
the pre-equilibrium phase of the medium prior to QGP 
formation. However, the strong chromo-electromagnetic 
field in this pre-cquilibrium stage of QGP may cause a 
Weibel Instability j3ll434| and influence the heavy quark 
thermalization process. Furthermore, the recent study 
[IHl indicates the radiation in the backward rapidity re- 
gion may have a larger impact on the heavy quark energy 
loss than expected before. 

In a follow-up work, we shall address some of these 
open questions, in particular how heavy quark ther- 
malization is affected by memory effects and the pre- 
equilibrium phase of the QGP evolution. 



ACKNOWLEDGMENTS 

We are grateful to the Open Science Grid for providing 
computing resources that were used in this work. We ac- 
knowledge many useful discussions with Berndt Miiller 
and Guangyou Qin and thank Christopher Coleman- 
Smith for supplying us with charm quark distributions 
from the PCM. This work was supported by U.S. de- 
partment of Energy grant DE-FG02-05ER41367. 



[1] M. Gyulassy and L. McLerran, Nucl. Phys. A750, 30 

(2005), nucl-th/0405013. 
[2] PHENIX Collaboration, S. Kelly, J. Phys. G G30, S1189 

(2004) , arXiv:nucl-ex/0403057. 

[3] STAR Collaboration, F. Laue, J. Phys. G G31, S27 

(2005) , arXiv:nucl-ex/0411007. 

[4] PHENIX, S. S. Adler et al, Phys. Rev. C72, 024901 

(2005) , arXiv:nucl-ex/0502009. 

[5] PHENIX, S. S. Adler et al, Phys. Rev. Lett. 96, 032301 

(2006) , arXiv:nucl-ex/0510047. 

[6] STAR, B. I. Abelev et al, Phys. Rev. Lett. 98, 192301 

(2007) , arXiv;nucl-ex/0607012. 

[7] PHENIX, A. Adare et al, Phys. Rev. Lett. 98, 172301 

(2007), nucl-ex/0611018. 
[8] PHENIX Collaboration, A. Adare et al, (2010), 

arXiv:1005.1627. 
[9] G. D. Moore and D. Teaney, Phys. Rev. C71, 064904 
(2005), hep-ph/0412346. 
[10] H. van Hees and R. Rapp, Phys. Rev. C71, 034907 

(2005), nucl-th/0412015. 
[11] Y. L. Dokshitzer and D. E. Kharzeev, Phys. Lett. B519, 

199 (2001), arXiv:hep-ph/0106202. 
[12] N. Armesto, A. Dainese, C. A. Salgado, and U. A. 
Wiedemann, Phys. Rev. D71, 054027 (2005), arXiv:hep- 
ph/0501225. 

[13] M. G. Mustafa and M. H. Thoma, Acta Phys. Hung. 

A22, 93 (2005), hep-ph/0311168. 
[14] M. G. Mustafa, Phys. Rev. C72, 014905 (2005), hep- 

ph/0412402. 



[15] M. Djordjevic and U. Heinz, Phys. Rev. C77, 024905 

(2008), arXiv;0705.3439. 
[16] Y. Akamatsu, T. Hatsuda, and T. Hirano, Phys. Rev. 

C79, 054907 (2009), arXiv:0809.1499. 
[17] P. Gossiaux and J. Aichelin, Phys. Rev. C78, 014904 

(2008) , arXiv:0802.2525. 

[18] B. Svetitsky, Phys. Rev. D37, 2484 (1988). 

[19] M. Golam Mustafa, D. Pal, and D. Kumar Srivastava, 

Phys. Rev. C57, 889 (1998), arXiv:nucl-th/9706001. 
[20] P. Gossiaux and J. Aichehn, J. Phys. G36, 064028 

(2009) , arXiv:0901.2462. 

[21] P. B. Gossiaux et al., (2011), arXiv:1102.1114. 

[22] H. van Hees, V. Greco, and R. Rapp, Phys. Rev. C73, 

034913 (2006), arXiv:nucl-th/0508055. 
[23] B. Combridge, Nucl. Phys. B151, 429 (1979). 
[24] C. Nonaka and S. A. Bass, Phys. Rev. C75, 014902 

(2007), nucl-th/0607018. 
[25] C. Shen, U. Heinz, P. Huovinen, and H. Song, Phys. Rev. 

C82, 054904 (2010), arXiv:1010.1856. 
[26] K. Geiger and B. MuUer, Nucl. Phys. B369, 600 (1992). 
[27] S. A. Bass, B. MuUer, and D. K. Srivastava, Phys. Lett. 

B551, 277 (2003), nucl-th/0207042. 
[28] M. Asakawa, S. A. Bass, and B. MuUer, Phys. Rev. Lett. 

96, 252301 (2006), hep-ph/0603092. 
[29] M. Asakawa, S. A. Bass, and B. MuUer, Prog. Theor. 

Phys. 116, 725 (2007), hep-ph/0608270. 
[30] M. Asakawa, S. A. Bass, and B. MuUer, J. Phys. G34, 

S839 (2007), arXiv:nucl-th/0702007. 
[31] E. S. Weibel, Phys. Rev. Lett. 2, 83 (1959). 



10 



[32] S. Mrowczynski, Phys. Lett. B314, 118 (1993). 125008 (2005), hep-ph/0408275. 

[33] P. Romatschke and M. Strickland, Phys. Rev. D69, [35] R. Abir, C. Greiner, M. Martinez, M. G. Mustafa, and 

065005 (2004), arXiv:hep-ph/0309093. J. Uphoff, (2011), arXiv;1109.5539. 

[34] P. Romatschke and M. Strickland, Phys. Rev. D71, 



